Equilibration rates and negative absolute temperatures for ultracold atoms in optical 

lattices 
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As highly tunable interacting systems, cold atoms in optical lattices are ideal to realize and ob- 
serve negative absolute temperatures, T < 0. We show theoretically that by reversing the confining 
potential, stable superfluid condensates at finite momentum and T < can be created with low 
entropy production for attractive bosons. They may serve as "smoking gun" signatures of equili- 
brated T < 0. For fermions, we analyze the time scales needed to equilibrate to T < 0. For moderate 
interactions, the equilibration time is proportional to the square of the radius of the cloud and grows 
with increasing interaction strengths as atoms and energy are transported by diffusive processes. 
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The concept of temperature is central in thermody- 
namics. For most systems, where the energy E has no 
upper bound, only positive temperatures, T > 0, are al- 
lowed in equilibrium. However, in any system with an 
upper bound in energy, T < are possible. In this case, 
states with higher energy are occupied more likely than 
states with lower E. Negative T and even phase transi- 
tions at T < have been realized in nuclear spin systems 
[H-Q: antiferromagnetically coupled nuclear spins show 
ferromagnetic order at T < where high-energy states 
are dominantly populated. Such an "inverted" popula- 
tion is also the basis of most lasers. 

While negative temperatures lead to many nonintu- 
itive results, all laws of thermodynamics can equally be 
applied Q . Figure HJa) shows schematically the entropy 
5 as a function of energy for a system with a maxi- 
mal and minimal energy. As 1/T = dS/dE, negative 
T arises whenever the entropy decreases as a function 
of energy. One consequence is that a Carnot engine, 
which operates between two reservoirs with temperatures 
T\ < and T 2 > 0, has an efficiency r\ larger than 1: 
V = W/Qi = 1 - T 2 /T 1 > 1, where W is the work done 
and Qi is the heat extracted from the (hotter) reservoir. 
Usually, rj < 1 as the entropy AS extracted from the hot 
reservoir has to be dumped into the second reservoir and 
the corresponding heat TAS is lost. If the first reservoir 
has, however, negative T, its energy is lowered when the 
entropy increases, see Fig. HJa): heat can be extracted 
from both systems simultaneously and therefore r] > 1. 

With quantum-optical methods it is possible to mani- 
pulate atoms in an unprecedented way. Atoms in optical 
lattices allows us not only to realize Hubbard models [B[ 
for bosons Q and fermions @, [H , but also to control all 
parameters with high precision. In the following, we will 
first show how simple manipulations of such systems can 
be used to create and detect unambiguously equilibrated 
quantum states at negative T, following partially a pro- 
posal of Mosk and then we discuss in the case of 
fermions the time scales needed for equilibration. 

Atoms in sufficiently deep optical lattices can be de- 



FIG. 1: (color online). (A) Entropy as a function of energy 
(schematically) for a system with an upper and lower bound of 
energy. For high energies, 1/T = dS/dE is negative. Arrows: 
for two reservoirs with T > and T < 0, respectively, one 
can remove reversibly energy from both reservoirs by reducing 
the entropy for T > and increasing it by the same amount 
for T < 0. (B) Schematic time-of-flight image of expanding 
bosons. For T > 0, repulsive bosons condense at momentum 
(with higher order peaks from Bragg reflections). For T < 0, 
condensates of attractive bosons form in the maxima of the 
kinetic energy for momenta (±ir/a,±n/a,±ir/a). 



scribed by Hubbard models [5] given by 
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for fermions and bosons, respectively. U is the local in- 
teraction, J the tunneling rate and a =f, \. the hypcrfine 
index for the fermions. Vo(t) represents a time-dependent 
parabolic potential with rii a = f^-f icr , Tii = bfb t . In 
Eqs. ([T][2]) higher bands are omitted. As discussed in de- 
tail by Mosk[9(, typical tunneling rates into such bands 
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are exponentially suppressed and negligible. 

Negative T in equilibrium is only possible for Hamilto- 
nians bounded from above. Thus, for the models ([IK]), 
Vb < and for bosons also U < is required. Never- 
theless, ultracold atoms have to be prepared with some 
Vb > initially, and we shall discuss how T < can be 
reached from such conditions (preparation of high-energy 
states in spin-systems is discussed in Ref. floj|). 

To get some intuition on negative T, note that the 
equilibrium density matrix, e ~ H / k BT ^ f or a Hamiltonian 
H at T < is identical to the density matrix for reversed 
temperature T = —T and a Hamiltonian H = —H. For 
the models (Q~|[2]) therefore considering T < is equiv- 
alent to T > with parameters — Vb, —U, and most 
importantly, — J. As for a cubic lattice with lattice con- 
stant a the sign of J can be absorbed into a shift of all 
momenta, k — > k + Q with Q = (ir/a,ir/a,Tr/a) [using 
11] that —Jcos(fca) = J cos(ka + it)}, the phase dia- 
gram of the negative- U Hubbard models for T < are 
identical to that of the positive-C/ Hubbard models if 
the momenta are shifted in all observables. Most dra- 
matically, bosons will therefore condense at momenta 
(±7r/a, ±7r/a, ±7r/a), i.e., in the maxima of the band 
structure, for U < 0,T < [ see Fig. QJb)]. Note 
that T < bosonic condensates are stable for attractive 
and unstable for repulsive U. For fermions superfluid- 
ity can also be reached for T < 0,U > but in this 
case the fermionic pairs condense at zero momentum as 
\52 k fk+Qff-k-QV We remark that it has 
12| that for [/ > J one can nevertheless 
induce a condensate at momentum Q even for fermions 
in a different nonequilibrium situation. To summarize, 
the momentum distribution of a Bose-Einstein conden- 
sate is probably the best way to detect T < due to 
the qualitative difference to the T > system, and we 
will quantitatively estimate under what conditions such 
a state can be reached. 

In principle, it is possible to reach T < without 
any entropy production by first adiabatically decreasing 
J until J = 0, then switching suddenly U — > —U and 
Vo — > — Vb, followed by an adiabatic increase of J un- 
til the desired value is reached. In practice this is not a 
good choice as the timescales for equilibration diverge for 
J — >• 0. Therefore one needs a faster scheme where J re- 
mains finite but as little entropy as possible is produced. 
As a reversal of Vq, U and T is formally equivalent to a 
sudden quench J — >• — J without inverting T (see above), 
the basic idea is to start from an initial state where J is 
finite but the kinetic energy is very small, i.e., a Mott or 
a band insulator. 

For bosons we propose the following. The system is (I) 
prepared in a Mott-insulating state. (II) J is switched 
off suddenly by increasing the intensity of the optical 
lattice. This freezes all density-density correlation func- 
tions. One waits (III) for a time t w during which the po- 
tential and the interaction are reversed slowly, Vq — > — Vq 
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and U — > — U. (IV) J is switched suddenly back to its 
initial value and (V) the system equilibrates. Finally, 
(VI) the negative trapping potential and/or the interac- 
tion strength U is weakened adiabatically with the goal 
to reach a superfluid condensate at T < 0. For fermions, 
parts of this scheme have been implemented in Ref. [Tljj , 
where, however, Vb was not reversed but set to 0. 

The waiting time t w allows us to reverse Vq and U 
slowly (manipulation of U requires "slow" changes of 
magnetic fields in some experimental setups), and more 
importantly, it leads to a complete dephasing of the ki- 
netic energy, E^i n , as has been argued in Ref. (ll|: Dur- 
ing tw each site in the lattice collects a different phase 
due to Vb- Effectively, therefore i?kin oc av- 
erages to zero for tw 3> 1/AV, where AV is the po- 
tential difference between neighboring sites. Note that 
for a Mott insulator E^m vanishes in the center of the 
trap but is sizable in the outer parts of the cloud where 
AV is large. Estimates using the parameters of Ref. @ 
show that waiting times i w > 10 ms are sufficient for a 
complete dephasing. 

For a quantitative estimate of the entropy generated 
during this sequence, one first needs to know the total 
energy E to t after step (V), i.e., before the final equilibra- 
tion stage. As the kinetic energy is zero after dephasing 
and the initial density distributions are fully preserved, 
we obtain E tot = —(E iai +E pot ), where E int and E pot are 
the initial interaction and potential energies, respectively. 
We obtain these variationally using the Gutzwiller wave 
function for bosons 13[ assuming a low-T initial state for 
a given strength of the confining potential characterized 
by the dimensionless compression VoV 2 / 3 /6J, where N 
is the number of atoms. This combination appears nat- 
urally |8| as the thermodynamic limit in a trap is given 
by N — > oo with V N 2 / 3 = const. 

Energy conservation allows us to determine the state 
after equilibration. The corresponding temperature will 
be negative and large. We therefore used a high- 



temperature expansion 14] to obtain the thermodynamic 
potential in the presence of the trap. We found that an 
expansion up to 4 th order in (3J gives accurate results. 
For simplicity, we locally approximated the system by 
a homogeneous one. Corrections to this so-called local- 
density approximation (li| vanish with 1/iV 1 / 3 and are 
therefore tiny for typical atom numbers, N ~ 10 5 . Fi- 
nally, the (negative) temperature was determined by set- 
ting (Hb) — Etot and then the entropy was computed 
from the high-T expansion. 

In Fig. [2] the entropy per boson is shown as a func- 
tion of the initial compression V V 2 / 3 of the cloud for 
several values of U. We also calculated the result with- 
out dephasing (dashed lines) , obtained in the limit tyj — > 
0, corresponding to the scheme originally proposed by 
Mosk Entropies and \T\ are about 40% higher in this 
case. 

The entropy per boson drops for increasing |Vb| as 
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FIG. 2: (color online). Entropy per boson as a function of the 
initial confinement Vo, for different initial interactions U after 
following the steps (I)-(V), see text. Dotted line: Condensate 
fraction at momenta (±7r/a, ±7r/a, ±7r/a) for U = 120 J af- 
ter a slow decrease of \U\ during step (VI). A slightly higher 
condensate fraction is obtained, if in step (VI) \Vo\ is reduced 
(not shown). Inset: inverse temperatures. Dashed lines: re- 
sults for U — 120 J using the protocol suggested by Mosk [t| , 
i.e., without dephasing of the kinetic energy. 



cloud where scattering is rare. In Ref. 
was also compared to experiments. 

We first consider an instantaneous quench, Voi —> 
Vof — — 0.05 Voi for U — 2J and parameters which reflect 
approximately the experimental conditions of Ref. [llj |. 
Temperature and density profiles (and parameters) are 
displayed in Fig. [3] Here T(r, t) is defined as the temper- 
ature of a homogeneous reference ensemble in equilibrium 
with the same energy and particle density. 
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larger and larger fractions of the cloud become Mott in- 
sulating and thus insensitive to sign changes of J. Most 
importantly, even for moderate compressions and mod- 
erate values of U the entropy per boson is well below the 
critical entropy of an ideal Bose gas in a d — 3 harmonic 
trap, S/N < sq — 3.6 ks- This implies that for realistic 
parameters, one can expect the formation of the negative- 
T Bose-Einstein condensate after the adiabatic expansion 
[step (VI)] characterized by the momentum distribution 
of Fig. [TJ3 and by a large condensate fraction. Figure [2] 
shows the condensate fraction assuming that in step (VI) 
\U\ has been reduced adiabatically to for fixed Vo < 0. 

We will now consider spin-1/2 fermions focusing on 
dimensions d = 2 to simplify numerics. With fermions 
one can reach T < using a time-independent U starting 
from a band insulating, rather than a Mott-insulating ini- 
tial state. We will focus our discussion on the time scales 
needed to reach T < 0. To determine the dynamics, we 
use a numerical solution of the Boltzmann equation in 
relaxation time approximation, 



dtf + Vk- V r / + F- V k / 



1 



r(n,e) 



(f-fo(n,e)), (3) 



where /(r, k, t) is the occupation probability in phase 
space and the force term F = — V r (Vbr 2 ) — UV r n(r) con- 
tains the external potential and interaction corrections on 
the Hartree level, n = n\ — ri\. = n{v,t) and e = e(r,t) 
are the local particle densities per spin and energy den- 
sity, respectively, and fo(n, e) is a Fermi function chosen 
such that energy and particle number are conserved. In 
Ref. ll| we have determined the scattering rate 1 /V(n, e) 
to reproduce the e- and n-dependent diffusion constant 
of the Hubbard model to order U 2 (obtained from an in- 
dependent calculation). Most importantly, l/r(n, e) van- 
ishes for low densities, l/r(n, e) oc n, in the tails of the 



FIG. 3: (color online). Cuts through the density n (lines) and 
inverse temperature f3 = l/kBT(r,t) (symbols) for quenches 
V 0i -)■ -0.05V 0i at U = 2J, V N = 60J, initial entropy 
S/N = 1.2fc s and N = 6000 particles per spin. 

Because of energy conservation and the upper bound 
on the kinetic energy, the atomic cloud cannot expand 
to arbitrary size but equilibrates for t — > oo. The initial 
T > becomes rapidly negative and slowly obtains a 
homogeneous value (which can be calculated from energy 
conservation in an independent way). It takes longer to 
equilibrate in the tails of the cloud, where scattering rates 
are small and strong Bloch oscillations occur. 

Two time scales determine the relaxation mainly. 
First, r(n, e) determines locally the equilibration rate ac- 
cording to Eq. ([3]) but local interactions do not change 
the local energy and particle density. Therefore a second 
time scale describes how long it takes to redistribute par- 
ticles and energy across the cloud. For not too weak in- 
teractions, see below, this transport will be diffusive and 
not convective as momentum is not conserved and can be 
efficiently transferred to the optical lattice by umklapp 
scattering for the relatively high \T\ considered here. The 
corresponding diffusion time is estimated as 



TD 



r 1 U Z N 
D 



J 3 



(4) 



where r is the radius of the cloud, D ~ v 2 t the diffusion 
constant, v ~ J a the typical velocity (for T ~ J), 1/r ~ 
nil 2 1 J and n{r/a) 2 ~ N. The inset of Fig. [4] displays 
1/T in the center of the cloud as a function of rrj/t. It 
shows that rjj becomes the relevant time scale for large 
N and not too small U. Furthermore, for t <C To, where 
the edges of the cloud do not yet play a role, there is 
a regime, where /3 ~ 1/t is observed as expected for 
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(energy) diffusion in d = 2. Figure [4] also shows the 
time t go {U) needed to reach J(3 = -0.469 (90% of the 
inverse temperature for t — > oo,U — > 0). For larger U 
one gets t 90 ~ t d oc U 2 but for small U, t 90 diverges 
due to the divergence of the local relaxation time, r oc 
l/U 2 . In all cases, relaxation is not very fast, tgo > 
300/ J. The fastest relaxation to equilibrium occurs for 
relatively weak interactions when U is a fraction of the 
full bandwidth 8 J. 




FIG. 4: (color online). Time required to reach /3q = 0.9 x 
/3 f (U -> 0) for TV = 6000. Inset: f3(r = 0,t) in the center 
of the trap as a function of rescaled 1/i for various U and 
system sizes (solid: TV = 3000, dashed: TV = 6000, dotted: 
TV = 12000). Dash-dotted line: J/3 = -0.469. 
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FIG. 5: (color online). The inverse temperature /3(r = 0,t) 
in the center of the cloud as a function of t/At for different 
opening times At, U = J and a Vo(t) shown schematically 
in the left inset {V of = -0.05V O i, V oi N = 60 J, TV = 6000). 
Right inset: Inverse temperatures /3(r = 0,t = oo) and /3(r = 
0, t = 2 At) as a function of At. 

In order to reach T < with small \T\, it is useful 
to decrease |Vo| slowly to reduce entropy generation and 
to see whether adiabatic conditions can be realized. We 
therefore use the protocol shown in the inset of Fig. [5] 
after a sudden quench at t = 0, Voi — > — Voi, —Vq is re- 
duced linearly, V (t) = -V oi + (V f + V 0l )t/At for t < At 
and Vo(t) = V f < for t > At. As shown in Fig. 
upon increasing At, considerably lower values of j3 < 
can be obtained and one approaches the adiabatic limit. 
Only due to the high entropy assumed for the initial state 



(S/N = 1.2 k B , implying Jp sa 1.47 for V oi -> \V of \ adi- 
abatically), |T| remains relatively high even for At — > oo 
where J/3 ~ —1.21. The overall entropy production for 
At — » oo is tiny, AS/N ps 0.12 kg, [c.f. the bosonic case, 
Fig. [5] as for the large initial T, kinetic energies were 
small. Note that even for At = 1600/ J deviations from 
the adiabatic behavior are considerable, which shows how 
difficult it is to reach truly adiabatic conditions. Never- 
theless, it is possible to reach T pb — 2 J within a time 
200/ J ~ 100 ms for typical parameters 
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In our opinion the observation of finite momentum su- 
perfluidity, Fig. QJ3, is probably the best "smoking gun" 
signature of T < in equilibrium. To reach it, it is, how- 
ever, necessary to switch the interaction U — > — U using 
a Feshbach resonance for bosons. We expect that the 
associated loss processes by three-particle scattering can 
efficiently be reduced in an optical lattice. An important 
issue are the time scales needed for local equilibration 
and - most importantly - redistribution of energy and 
particles across the system. For fermions we find that re- 
laxation is most efficient for relatively weak interactions. 
More generally, the long equilibration times arising from 
the necessity to redistribute energy and particles should 
be important for all equilibration processes and quenches 
in inhomogeneous systems both for positive and negative 
T. Here we expect that the equilibration properties of 
low- 1 T | bosons differ qualitatively from high-|T| fermions 
due to the suppression of umklapp scattering for bosons 
and due to the onset of superfluidity. 
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E. Demler, and U. Schneider, and financial support by 
the SFB 608 and SFB/TR 12 of the DFG and the 
Studicnstiftung des deutschen Volkes (S.M.). 
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